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Abstract. We present an algorithm for the symbolic and numerical compu- 
tation of the degrees of the Chern-Schwartz-MacPherson classes of a closed 
subvariety of projective space PJ5. As the degree of the top Chern-Schwartz- 
MacPherson class is the topological Euler characteristic, this also yields a 
method to compute the topological Euler characteristic of projective varieties. 
The method is based on Aluffi's symbolic algorithm to compute degrees of 
Chern-Schwartz-MacPherson classes, a symbolic method to compute degrees 
of Segre classes, and the regenerative cascade by Hauenstein, Sommese and 
Wampler. The new algorithm complements the existing algorithms. We also 
give an example for using a theorem by Huh to compute an invariant from 
algebraic statistics, the maximum likelihood degree of an implicit model. 



1. Introduction 

The topological Euler characteristic is a widely studied invariant in many areas 
of mathematics. In this article, we present an algorithm for the symbolic and 
numerical computation of the topological Euler characteristic of a possibly singular 
closed subvariety of Pg. More generally, we present a method to compute the 
degrees of the Chern-Schwartz-MacPherson classes of such a variety. As the degree 
of the top Chern-Schwartz-MacPherson class is the topological Euler characteristic, 
this yields a method to compute the topological Euler characteristic for projective 
varieties. 

The method we describe is based on several different known methods. According 
to Aluffi's article [2J, degrees of the Chern-Schwartz-MacPherson classes can be 
computed from so-called shadows of graphs, also called projective degrees of a 
map. Using other results from [2], we prove in this article that the shadows of 
graphs can be computed from the degrees of Segre classes, which in turn can be 
computed from so-called residuals, according to [5]. The residuals can be computed 
symbolically, but are also a by-product of the regenerative cascade, a method for 
computing the numerical irreducible decomposition of a variety developed in [T3] . 

The existing methods to compute the topological Euler characteristic of varieties 
are, to the author's knowledge, all symbolic. The algorithm developed in [2] was im- 
plemented by Alum in a Macaulay2 package called CSM. Furthermore, Macaulay2 
provides the command euler which computes the Euler characteristic using 
the computation of Hodge numbers. However, the latter only works for nonsin- 
gular varieties. The algorithm described here is implemented in version 0.2 of the 
Macaulay2 package CharacteristicClasses. More details about the implementation 
can be found in |15j . 

We give a short outline of the content of this article. In Section 2, the definitions 
of Chern-Schwartz-MacPherson classes, Segre classes and their degrees are recalled. 
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In Section 3, we describe how to compute degrees of Chern-Schwartz-MacPherson 
classes from the degrees of Segre classes, and in Section 4 how to compute degrees 
of Segre classes numerically using the regenerative cascade. In Section 5, the run- 
ning times of our implementation are compared to existing implementations. We 
conclude with an example from algebraic statistics in Section 6. 

2. Segre classes, Chern-Schwartz-MacPherson classes and their 

DEGREES 

We work over the field of complex numbers. In this and the following section we 
will use the language of schemes instead of varieties. Let X be a possibly singular 
closed subscheme of P™, embedded by i : X c — > P™, and let k be the dimension of X. 
We will recall the definitions of the Chern-Schwartz-MacPherson classes and Segre 
classes of X as well as their degrees. The standard reference for intersection the- 
ory, including Chow groups and Segre classes is [TO]- Chern-Schwartz-MacPherson 
classes are described in several articles by Alum, e.g. [2J. 

Characteristic classes, such as Chern-Schwartz-MacPherson and Segre classes, of 
X are elements of the Chow group A„{X) = ®j =0 Aj(X) of X, i.e., they are cycles 
modulo rational equivalence. However, the Chow group may be hard to compute 
in general. Hence computational methods for characteristic classes focus on com- 
puting the degrees of the respective classes rather than the classes themselves. Let 
a = YH=o a i[Vi] be a d-dimensional cycle in Ad(X), then its degree is defined as 
deg(a) = X)i=o a * deg(T^), where deg(Vi) is the degree of Vt considered as a subva- 
riety of P n . Equivalently, one can consider the pushforward i*(a) of a to the Chow 
group A*(P n ) of P". Let H be the hyperplane class, then A*(P") = Z[H]/(H n+1 ) 
and i*(a) = H n ~ d Y^i=o a «deg(Vi). Observe that the definition of degree does not 
agree with the usual definition of the degree of cycle classes denoted by J . 

The total Chern class of a nonsingular variety is defined to be the total Chern 
class of its tangent bundle. There are several generalizations of this concept to possi- 
bly singular schemes, including Chern-Schwartz-MacPherson classes. They were de- 
fined independently by MacPherson [TB] , proving a conjecture of Grothendieck, and 
Schwartz |18) . and were shown to agree in [5]. Chern-Schwartz-MacPherson classes 
enjoy nice functorial properties, which we recall by summarizing sections 2.2 and 2.3 
of [2J. Let S be a proper scheme. The Chern-Schwartz-MacPherson class of a closed 
subscheme X of S is an element csm(X) in A*(X) such that csm(X) = c(X) (~l [X], 
the total Chern class of X, for nonsingular schemes X. The i-th Chern-Schwartz- 
MacPherson class (csM)i(X) 6 Ak-i(X) is then the codimension i part of the total 
Chern-Schwartz-MacPherson class csm(X). Chern-Schwartz-MacPherson classes 
extend to constructible functions by csm(X)vcs m v^-v) = XVcS m vcsM(V). This 
actually gives a natural transformation C ~+ A from the functor of constructible 
functions C to the Chow group functor A. The functor C of constructible functions 
maps a scheme to the abelian group of constructible functions on it. A morphism 
/ : S — > T of schemes is mapped to the morphism C(f) of abelian groups by 
C(/)(ly)(y) = x(/ _1 (y) n V"), for V C S a subscheme and y £ T a closed point. 
Here x 1S t ne topological Euler characteristic and hence x{f~ l {y) H V) the topo- 
logical Euler characteristic of the fiber of the point y. As a special case of the 
functorial properties, we get that Chern-Schwartz-MacPherson classes compute the 
Euler characteristic of the support of schemes. Let n: X — > point. Then 

k*c sm (A) = c S m{C(k)1x) = csM(x(A r cd)lpoint) = x(A rcd ) [point] , 

hence J csm(A) = x(A re( j), where J as usual denotes the degree of the top class. 
Observe that constructible functions follow laws of exclusion-inclusion, e.g., 

lxi nx 2 — lxi + lx 2 — lx 1 ux 2 < 
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It follows that Chern-Schwartz-MacPherson classes follow similar laws of exclusion- 
inclusion, e.g., 

csm(^i n X 2 ) = c SM (Xi) + csm(X 2 ) - c S m{Xi U X 2 ). 

Finally, we recall the definition of Segre classes of closed subschemes of P™ . The 
i-th Segre class Si(X,¥ n ) of X embedded in P™ is defined to be the i-th Segre 
class of the normal cone CxP" of X in P". Instead of using the definition of 
Segre classes of cones, we use a shortcut, allowing us to define the Segre classes 
Si(X,F n ) G Ak-i for < i < k directly. Let Ex be the exceptional divisor of 
the blow-up BlxP" of P™ along A, and let r\: E x — > X be the projection. Then 
Si(X,P n ) := ( — l) p 77* (£7^) for p = n + % — k, where E x denotes the p-th self 
intersection of the exceptional divisor Ex- The total Segre class s(A, P") of X in 
P" is then defined as s(X, P") = 1 + si(X, P n ) + . . . + s k (X, P"). 

3. Computing degrees of Chern-Schwartz-MacPherson classes using 

degrees of segre classes 

In this section we combine results from [2] with a theorem proved here to give 
an algorithm for the computation of the degrees of Chern-Schwartz-MacPherson 
classes, provided an algorithm for the computation of degrees of Segre classes. 
As in the previous section, we consider a fc-dimensional closed subscheme X of 
P™, embedded by i: X P". Recall that computing the degrees of the Chern- 
Schwartz-MacPherson classes is equivalent to computing the pushforward ?*csm(A) 
of the total Chern-Schwartz-MacPherson class to the Chow group A*(F n ) of P™. 

Observe first that it suffices to find an algorithm computing the Chern-Schwartz- 
MacPherson classes of hypersurfaces. The Chern-Schwartz-MacPherson classes of 
lower-dimensional schemes can then be computed using the exclusion-inclusion prin- 
ciple described in the previous section: 

i*csM(Xi n X 2 ) = i*csm(Xi) + i*c SM (X 2 ) - i*c S M(Ai U X 2 ). 

Due to a result from [2], the Chern-Schwartz-MacPherson classes of a hypersurface 
can be computed from the so-called shadow of the graph of the singular locus of 
the hypersurface. Reversing another result from [5J, we prove here that this shadow 
can be computed from the degrees of the Segre classes of the singular locus. We 
shortly review the definitions and results from [5J. 

Let X be a hypersurface in P™, given by a homogeneous polynomial / in the 
polynomial ring C[xq, . . . ,x n ]. The singular locus of X is given by the ideal 
{■g~-, ■ ■ ■ , ■§§-)■ The blow-up of P™ along the singular locus of A" is a closed 
subscheme of P" x P™ which also can be seen as the closure of the graph of 
(-§f : . . . : : P n ^ ¥ n . Denote its class in the Chow ring of P™ x P" by 
r. The structure theorem for the Chow group of projective bundles (Theorem 3.3 
in [TU]) assigns to T a cycle in the Chow group of the first factor P™. Aluffi calls this 
cycle the shadow of T, and denotes it by G = Y^i=o 9i^ 1 - The numbers go, . . . ,g n 
are also called the projective degrees of the map (J^- : . . . : ■§£-), see chapter 19 of 
[12]. According to Lemma 4.2 in [3], the graph of T can be computed by intersecting 
r with the pullback of the hyperplane class from the second factor and projecting 
to the first factor. The importance of the shadow of the graph of the singular locus 
is that it can be used to compute the degrees of the Chern-Schwartz-MacPherson 
and Segre classes of A, according to the following propositions, Theorem 2.1 and 
Proposition 3.1 from [2J. 
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Proposition 1 (Aluffi). Using the notation above, the pushforward of the Chern- 
Schwartz-MacPherson class of X is given by 

n 

^c SM (X) = (1+H) n+1 -j2 gj (-H)S(l + H)»-i. 

3=0 

Proposition 2 (Aluffi). Let Y be a subscheme o/P™ given by an ideal with gen- 
erators all of the same degree r, and let G be the shadow of the graph ofY. Then 
the pushforward of the total Segre class of Y is given by 

i*s{Y, P") = 1 - ciOirH))- 1 n (G ® 0(rH)). 

According to Definition 2 in pQ, G ® 0{rH) = £™ =0 c(S(rH)y • 

We prove a theorem reversing Proposition [5J i.e., a method to compute the 

shadow of the graph of the singular locus from the degrees of its Segre classes. 

Together with Alum's results and a method to compute the degrees of the Segre 

classes of a projective scheme, this yields a method to compute the degrees of the 

Chern-Schwartz-MacPherson classes of a projective scheme. 

Theorem. Let X be a hypersurface in P™. Let so,...,s^ be integers related to 
the degrees of the Segre classes of the singular locus of X by Sq = 1, Si = . . . = 
Sn-fe+l = and s~i = — deg Sj_( n _fc)(sing(X),P n ) for i = n — k, . . . ,n, where 
k = dim(sing(A)). Let the shadow of the graph of sing(A) be denoted by G — 
go + . . . + g n H n . Then the integers gi are related to the integers Si by the following 
triangular linear equation system with ones on the diagonal: 



E K-r^ 



3=0 

where i = 0, . . . ,n and r is the degree of the generators of the ideal of sing(A). 
Proof. Denote Y — sing(A). According to Proposition^ it holds that 

Us(Y, P) = 1 - c(C(riJ))- 1 n (G (8) 0{rH)). 
Recall that G ® O(rH) is defined to be YT l= o c(o(rH)Y ■ The identity can be rewrit- 



ten as 

i = '4 = ^ ' 1 = V ' 

We develop the right-hand-side of this identity: 

i=0 v ; i=0 \j=o I i=0 i=0 v 7 

where the last step uses the following identity for formal sums in a variable x: 

i+l 



I 



X 



5> =E 

j=0 / 1=0 

for any positive integer i. We continue rewriting the right-hand-side of the equation: 

e^eh V^ = 

i=0 i=0 v ' 
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Summarizing, we have proved that 



Comparing the coefficients of the powers of H on both sides yields the equations 
to be shown. □ 

Corollary. Let X be a hypersurface in P™ . Given the degrees of the Segre classes 
of the singular locus of X and the integers So,...,s^ related to them as in the 
theorem before, the degrees of the Chern-Schwartz-MacPherson classes of X can be 
computed in the following way. Solve the linear equation system 



3=0 

for the go, . . . ,g n . Then compute the Chern-Schwartz-MacPherson classes by 

deg (c SM ) p (X) = ("t 1> L )'(" •)• 
where i = n — dim(X) + p and p = 0, . . . , dim(X) . 

Proof. The first part is the theorem above. The second part is a reformulation of 
Proposition [H which states that 

n 

i*c SM (X) = (1 + H) n+1 - 9l{-H)\l + H) n ~ l . 

1=0 

We rewrite the second term in the right-hand-side of the identity: 



n — I 



Hence 



^ gi (-H) l (l + HT- 1 = ]T gi(-H) 1 JJ(" ")H* 

n n—l / 7\ n i ✓ 

;=o k=0 V / i=0 j=0 v ■' 



z*c S m(^) = 




Comparing the coefficients of the powers of H on both sides yields the identity to 
be shown. □ 

Recall that it suffices to compute the degrees of the Chern-Schwartz-MacPherson 
classes for hypersurfaces due to the inclusion-exclusion principle. We hence have 
found an algorithm to compute the degrees of the Chern-Schwartz-MacPherson 
classes and hence the topological Euler characteristic of any closed subscheme of 
P", provided an algorithm for the computation of Segre classes. 

Observe that the computation of the topological Euler characteristic of projective 
schemes also allows to compute the Euler characteristic of affine schemes. One just 
computes the difference between the Euler characteristic of the homogenization and 
the subscheme at infinity. 

Example. We compute the degrees of the Chern-Schwartz-MacPherson classes of 
the nodal plane cubic given by the ideal (x 3 + x 2 z ~ y 2 z) C C[x, y, z\. The singular 
locus of the curve is the point P = [0 : : 1] and the degree of its Segre class 
so(P, P 2 ) is the degree of the point, deg(so(-P, P 2 )) = 1. Hence §q = 1, s[ = and 
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S2 = — 1 • The shadow of the graph of the singular locus can thus be computed via 
the following linear equation system 

1 = s~ = Q (-2)°-% = go 

= Si = Q (-2) 1 -% + Q (-l) 1 " 1 ^ = (-2) ff0 + 5i 

-1 = *~2 - Q (-2) 2 -% + (J) (-2) 2 - 1 .g 1 + Q (-2) 2 ' 2 52 = 4 5o - 4 fll + ,g 2 . 

Hence go = 1, <?i = 2 and 172 = 3. Now we can compute the degrees of the Chern- 
Schwartz-MacPherson classes of the singular plane cubic C. 



deg(c SM )o(C) 



r)-*-<:S)-«<-'>'(:::)- 

deg(c SM )i(C) = ( 2 2 ') - 9o(-l)° Q I 0) - ^i(-l) 1 (2 1 1) " 9 ^ 



2-2 
2-2 



1. 



This confirms that the nodal plane cubic has topological Euler characteristic 1. 
The pushforward to the Chow ring of P 2 of the Chern-Schwartz-MacPherson class 
is «*csm(C) = 2>H + H 2 , where H is the hyperplane class. 



4. Computing degrees of Segre classes using the regenerative 

CASCADE 

In this section we describe how to compute degrees of Segre classes numerically. 
The main ingredience of the method is the algorithm in [5] , a symbolic algorithm for 
the computation of the degrees of Segre classes. In more detail, the method in [5] 
reduces computation of degrees of Segre classes to the computation of the degrees 
of so-called residuals, which can be computed symbolically. In the following, we 
describe how the regenerative cascade from [13] can be used to compute the resid- 
uals, which was already claimed in The regenerative cascade is implemented 
in Bertini [¥] , a software for the numerical solution of polynomial equation systems 
via homotopy continuation. One of the advantages of homotopy continuation is 
that the algorithms are parallelizable. 

We start by recalling the Segre class algorithm from [9]. Let X be a closed k- 
dimensional subvariety of P™ given by an ideal /. Let m be the maximal degree of 
the generators. Pick n random degree m elements /1, . . . , /„ of the ideal /, corre- 
sponding to n hypersurfaces containing X . According to a Bertini type theorem, d 
of these hypersurfaces intersect in the scheme X and a residual scheme Rd of codi- 
mension d, where d = n — k, . . . ,n. Furthermore, one obtains the following relations 
between the degrees of the Segre classes of X and the degrees of the residuals 

deg( Sp (X,P")) = m d - deg(i? rf ) - £ _ .jm^Meg ( Sl (X,P n )) , 

i=o \P l ' 

for p = d — (n — k). One observes that these relations form an upper-triangular 
linear equation system with ones on the diagonal. Hence the degrees of the Segre 
classes can be computed easily once the degrees of the residuals are known. 

We continue with describing the regenerative cascade developed by Hauenstein, 
Sommese and Wampler in 13J. It is a method to compute the so-called numer- 
ical irreducible decomposition of a variety, i.e., to compute (and represent in a 
specified way) the components of the solution set of a polynomial equation system. 
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As we will see, the degrees of the residuals are a by-product of these computa- 
tions. Let a polynomial equation system be given by homogeneous polynomials 
/i> • • • > fn £ C[xo,. .., x n ] of the same degree m. The aim of the computations is a 
set of points containing so-called witness sets for each component of the numerical 
irreducible decomposition of V(fi, . . . ,/„) C P™. A witness set for a component 
is the intersection of the component with a general linear space of complementary 
dimension, hence a number of generic points on the component. In general, the 
polynomial equation system need not have as many equations as the dimension of 
the ambient space and the equations can have different degrees. However, this spe- 
cial case is sufficient for our purposes. Let L\, . . . , L n be random linear functions in 
C[x , ...,x n ] and let T k = . . . , fk,L k +i, • • • , L n ] T for k = 1, . . . , n. A point in 
V(J-fe) is called a non-solution if it does not lie in V(fi, . . . , /„). The non-solutions 
are nonsingular and isolated. The number of non-solutions is a by-product of the 
regenerative cascade computations. Geometrically, all solutions of Tk are obtained 
by intersecting the variety V(/i, . . . , fk) with a general linear space of codimension 
k. If the polynomials /i, . . . , /„ are chosen in the same way as for the Segre class 
algorithm above, the non-solutions are the intersection of the residual with the lin- 
ear space. As the residual is pure-dimensional, the number of non-solutions is the 
degree of the residual. The by-product of the regenerative cascade, the number of 
non-solutions of J-fc , is hence exactly the information needed to compute the degrees 
of the Segre classes of X . 

Example. The twisted cubic Ct w is given by the ideal (xz — y 2 , yw — z 2 , xw — yz) C 
C[ar, y, z, w]. We randomly choose three polynomials 51,32,53 of degree r = 2 in 
the ideal. By using the implementation of the regenerative cascade in Bertini [4], 
we compute that the number of non-solutions of [51, 52, L^\ T is deg(i?2) = 1, and 
the number of non-solutions of [5i,52,53] T is deg(i?3) = 0. This is consistent with 
the fact that two general degree 2 elements of the ideal of the twisted cubic cut out 
the twisted cubic and a line, and that three general degree 2 elements cut out only 
the twisted cubic. Hence the degrees of the Segre classes of the twisted cubic can 
be computed in the following way: 

deg(s (C tw ,P 3 )) = m d - deg(i? 2 ) = 2 2 - 1 = 3, 

deg( Sl (C tw , P 3 )) = m d - deg(i? 3 ) - (f\m ■ deg(s ) - 2 3 - - 3 ■ 2 • 3 = -10. 

5. Implementation and comparison to other implementations 

As said above, the algorithms described in this article are implemented in the 
Macaulay2 [11] package CharacteristicClasses. The implementation is described in 
a more detailed way in [15]. 

We compare the running times of our implementation to Aluffi 's implementation 
CSM, December 2011 version, and the routine euler from Macaulay2. Observe 
that the latter only works for nonsingular varieties and takes a projective variety 
as input. The command CSMclass from our implementation uses the algorithm 
described here. The command CSM from Aluffi's implementation with the same 
name uses the algorithm described in [2], and the command euler computes the 
topological Euler characteristic from the Hodge numbers. We use examples similar 
to the ones in [S], and a two core processor with 1.40GHz and 4MB RAM. The 
results are summarized in Table [T] 

The equations defining the different examples are found on the author's webpage, 
www. math. su. se/^jost/examplesCSM.m2. All computations were done over the 
rational numbers. The ideal of the smooth surface is generated by the 2-by-2 minors 
of a 2-by-3 matrix of random linear forms. 
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Table 1. Comparison of run times. 



Input 


CSMclass 


CSMclass 


CSM 


culcr 




symbolic 


numeric 






twisted cubic 


< Is 


47s 


2s 


< Is 


smooth surface in P 4 defined by 


34s 


2131s 


14777s 


88s 


minors 










Segre embedding of P 1 x P 2 in P 5 


10s 


285s 


Is 


< Is 



As Table [T] shows, the different symbolic implementations complement each 
other. As the run time of both CSMclass and CSM is exponential in the num- 
ber of generators, the command euler may be the best choice for examples with 
many generators. On the other hand, euler only works for smooth varieties. The 
numerical implementation is actually slower than the symbolic implementation for 
small examples such as shown in Table [TJ However, for very large examples the 
symbolic methods may not terminate due to insufficient memory. The numerical 
implementation does not have this problem. Furthermore, the numerical algorithm 
is parallelizable as the main step, the regenerative cascade, is parallelizable. 

The complexity of computing the topological Euler characteristic has been stud- 
ied in [BJ and been found to be complete in the class FP^ Pc . Roughly speaking, 
functions in #Pc are polynomial time counting functions and functions in FP^ Pc 
are functions computable in polynomial time using oracle calls to functions in #Pc- 
For the precise definition of the class, we refer to, e.g., [BJ. 

6. Example from algebraic statistics 

We illustrate the algorithm for the numerical computation of the Euler charac- 
teristic by an example from algebraic statistics which uses the topological Euler 
characteristic: the maximum likelihood degree, defined in [7J. We compute the 
maximum likelihood degree of the statistical model described in Example 2.2.2 in 
[8] using a theorem from |14) together with our implementation of the algorithm 
presented here. 

To introduce the maximum likelihood degree properly is beyond the scope of 
this article. The reader interested in details may consult [5] or [17] , 

In algebraic statistics, a statistical model for discrete data is given by a subvariety 
X of P™ = Pg. In this context, we focus on implicit models, i.e., subvarieties given 
by equations. Denote by poi ■ • ■ ,Pn the homogeneous coordinates of P™, where pi 
represents the probability of the i-th event. Of course, the statistically interesting 
case is that all coordinates are real and positive. A data vector is given by natural 
numbers uq, ■ ■ ■ ,u n , where Ui is the number of times the i-th event was observed. 
The problem of maximum likelihood estimation is to find the point in the model 
which best explains the data vector u. This is done by maximizing the likelihood 
function 

L(p Q ,...,p n )=p%° -...-p^. 

One usually assumes that none of the probabilities vanishes, and that they sum up 
to 1. Hence one usually maximizes the likelihood function on the open subvariety 
U given by 

U = {x £ X\p ■ ...-p n (pi + ...+p n )¥ : 0}- 
The likelihood function can in general have several critical points. Hence methods 
finding only local maxima such as the Newton method and its derivates may run 
into problems. It thus makes sense to define the maximum likelihood degree of a 
model to be the number of critical points for a general choice of the data vector 
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u. The maximum likelihood degree was introduced in [7J. Since then considerable 
amount of research has been done computing the maximum likelihood degree of 
certain classes of models. 

In |14j , Huh generalizes a conjecture of Varchenko on the maximum likelihood 
degree of hyperplane arrangements to smooth very affine varieties. A very afhne 
variety is one that can be embedded into an algebraic torus. The following theorem 
is a special case of Theorem 1 in [T3] . 

Theorem (Huh). If the exponents it, are sufficiently general, then the number of 
critical points of the likelihood function is equal to the signed Euler characteristic 

Example. We use this theorem on the random censoring model with two events 
described in Example 2.2.2 in |S]. It is given by the ideal / = (2p pip2+plp2+PiP2~ 
p\p\2 +P1P2P12) in the polynomial ring C\po,pi,P2,Pi2\- With our implementation 
of the algorithm described here, we compute the Euler characteristic of V(I) to be 
5. As we are interested in the Euler characteristic of U = V(I) \ V (poPiP2Pi2(po + 
Pi +P2 + P12)), we compute the Euler characteristic of V {poPiP2Pi2(po +Pi +P2 + 
P12)) H V(I), which is 2. So the Euler characteristic of U is 5 — 2 = 3, and we get 
that the maximum likelihood degree of the model is 3. This agrees with the result 
in [5] . The computations for this example are presented in a more concrete way in 
Example 3 in [15] . 
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